Kinetics of the A + B — > reaction with mass-dependent fragmentation 
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We investigate the kinetics of uniformly driven A + B — ► reaction with mass-dependent frag- 
mentation in one dimension. In this model, the fragmented mass m of a site with mass m is given as 
m = Ui , and it is driven to the one direction. When opposite species masses occupy the same site, 
mass reaction takes place instantaneously. Since the fragmented mass m of A < 1 is less than mass 
Hi of a site, the exponent A controls the attractive interaction between particles at the same site. 
The A = case corresponds to hard-core (HC) particle system. With equal initial densities of both 
species, we numerically confirm that the scaling behaviors of density and lengths except the domain 
length I are the same as that of the uniformly driven HC particle system. The scaling behavior of 
t is t ~ t 2 ^ 3 . The kinetics of the reaction is independent of A as long as A < 1. The A-independent 
kinetics results from the A-independent collective motions of single species domains. 

PACS numbers: 05.70.Ln,05.40.-a 



I. INTRODUCTION 

The irreversible two species reaction A + B — > has 
been widely and intensively investigated among diffusion- 
limited reactions because of its rich kinetics and wide ap- 
plications to various phenomena such as physics, chem- 
istry and biology [J S, I, 1 0, i . The reaction instanta- 
neously takes place with a rate k when two particles of 
opposite species encounter on the same site. For homoge- 
neous initial distributions with equal densities of A and 
B, pa(0) = Pb(0), the density p(t) algebraically decays 
in time t. In higher dimensions than the upper critical 
dimensions (d > d c ), the kinetics follows the mean-field 
rate equation and p decays as p ~ t^ 1 in time t. How- 
ever for d < dc, the kinetics is fluctuation-dominated 
and depends on the motion and the mutual statistics of 

particles @, 1, @, 0, S ©, EE ED, El HI 0, HE HE Ezl . 

For isotropic diffusions, p(t) scales as p(t) ~ t~ d / 4 with 
d c = 4 [E IE 0, H, Hi ■ However, when the diffusion of two 
species is anisotropic, the kinetics is completely changed. 
With the global relative drift of one species to the other, 
p(t) scales as p(t) ~ £-( d+1 )/ 4 for d < 3 @. When only 
boundary particles of domains are relatively driven, the 
kinetics continuously changes according to the form of 
the drift velocity [lfj. For all cases mentioned above, 
the hard-core (HC) constraint between identical particles 
is irrelevant to the asymptotic kinetics. However, when 
both species are uniformly driven to the same direction, 
the HC constraint completely changes the kinetics. In 
one dimension, Janowsky reported that p scales as i -1 / 3 
for HC particles instead of t~ x l A [TtJ] . Since one expect 
that the motion of particles under the Galilean transfor- 
mation is isotropic, the t~ d ^ 4 decay law of the isotropic 
diffusion is also expected in uniformly driven systems. 
Furthermore, for the single species reaction kA — ► 0, the 
anisotropy of diffusion is no effect on the kinetics of HC 
particles [19j. Hence the result of Janowsky is rather 
surprising. For the uniformly driven A + B — > reaction 
without HC exclusion, the uniform drift also does not 
change the kinetics [2fJ. Hence the interplay of the uni- 



form drift and HC constraint is crucial for the anomalous 
i -1 / 3 decay law. 

As an attempt to understand the origin of the anoma- 
lous decay law i -1 / 3 , Ispolatov et al suggested a scaling 
argument provided the motion of particles in a single 
species domain follows the Burgers equation [ijj]. From 
the drift velocity of a single interface between A and B 
domains, they found the rate of the change of the do- 
main length I, defined as the distance between the first 
particles of adjacent opposite species domains, scales to 
the density p as dt/dt ~ p. For random initial conditions 
of pa{0) = Pb{0), P scales as p ~ \j\fl 6} so I scales 
as £ ~ i 2 / 3 instead of i 1 / 2 of the isotropic case. Hence, 
the anomalous scaling of the domain length £ may re- 
sult in the anomalous density decay. In d dimensions, I 
scales as I — £( 5 - d )/ 6 . So p decays as p ~ £-( d+1 )/ 6 for 
d < 2, t- d ^ for 2 < d < d c (= 4) 0. In d > 2 dimen- 
sions, the drift is irrelevant. In addition to the domain 
length £, the inter-particle distance (Iaa) and the distance 
between two adjacent particles of opposite species(£4R) 
also characterize the spacial organization of particles [7[ . 
Janowsky numerically showed that £aa and £ab scale as 
Iaa ~ t 1 ^ 3 and £ab ~ t 3 ^ 8 respectively [T^]. Intriguingly, 
the scaling of £ A b is not changed by the drift. Using the 
scaling of £ab, it was suggested that the domain length 
I scales as i 7 / 12 [13j |. The extensive simulation results of 
Ref. [IH supported the result of Janowsky, but it con- 
flicts with the prediction of Ref. [l2j . 

Recently there have been some attempts to describe 
systems of HC particles in field-theoretic formalism 
[lE HI HI- S.-C. Park et al systematically derived the 
i -1 / 3 decay law of p. The scaling of p is well under- 
stood, while the scaling behaviors of various lengths are 
not convincingly conclusive. Since the asymptotic scal- 
ing regime is extremely slowly approached due to strong 
corrections [HI, [3| , it is hard to observed the true scaling 
behaviors directly via simulations. Hence it is desired to 
study another model which exhibits the same scaling be- 
havior, but reaches the true scaling regime within moder- 
ate simulation time. On the other hand. HC interaction 
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FIG. 1: Motions of particles driven to the right, (a) Motions 
of HC particles. Only boundary particles (gray boxes) can 
hop to the right, (b) Motions of bosonic particles with A = 0. 
A randomly selected particle (gray box) hops to the right. 



can be generalized by allowing multiple occupation with 
a certain attractive interaction between particles on the 
same site. The on-site attractive interaction would give 
more general understanding about the scaling behavior 
HC particles exhibit. 

For this aim, in this paper, we introduce a one di- 
mensional two species bosonic reaction model with mass- 
dependent fragmentation. In this model, the HC con- 
straint is replaced by mass-dependent fragmentation of 
bosonic particles. For simplicity, we consider the driven 
case to the right direction only. If the randomly selected 
site i is occupied by n, identical particles or mass n, , the 
mass m — n\ moves to the right nearest neighbor with 
unit rate. The fragmentation of mass to corresponds to 
the moves of randomly selected m particles among m 
particles. When the opposite species of mass n$+i oc- 
cupies the site i + 1, the reaction — > n,j+i — m 
occurs instantaneously. For to > ni+j., to — rij+i mass 
occupies the site i + 1. For A — 1, the whole mass rii 
moves to the right so A = 1 corresponds to the uniformly 
driven bosonic model. For A < 1, the fragmented mass is 
less than m, which can be regarded as the consequence 
of the attractive interaction between particles occupying 
the same site. Our model generalize the HC interaction 
and enables to investigate how the kinetics of the reaction 
depends on the interaction strength. 

Especially, for A = 0, only one particle or unit mass 
hops to the right with unit rate. For HC particles, since 
only the right boundary particle of a train of particles can 
hop with unit rate, the A = case corresponds to HC in- 
teraction (Fig. [T|) . In some bosonic lattice gas model in 
which at each site, only the top particle on the pile of 
particles hops to the bottom of the pile of the nearest 
neighboring site. In this model, the order of particles is 
conserved so we call this model the ordered bosonic lat- 
tice gas (BLG). It was shown for the drift case that the 
motion of particles follows the Burgers equation as in HC 
particle systems [23[. Since particles on a site are iden- 
tical, the ordered BLG describes the motion of particles 
inside single species domains of the A = case in our 
model. For this reason, we regard HC interaction as the 
on-site attractive interaction between bosonic particles 
in our model. Hence the kinetics of A = is expected 
to be the same as that of HC particles. For < A < 1, 



the number of hopping particles from a site (to) is larger 
than one, but still smaller than m. Hence the interaction 
is weaker than HC constraint. However, the interaction 
of A < is stronger than HC interaction due to m < 1. 

We study the kinetics of the reaction A + B — > with 
mass-dependent fragmentation rate D{m) = m x for ran- 
dom initial distributions of equal densities of A and B. 
With uniform drift to the one direction of both species, 
we find for several values of A that the scalin g of p, Iaa 
and Iab coincide with the results of Ref . 13l Il4j . How- 
ever the domain length t scales as I ~ t 2/ ^of Ref. [l2j 
rather than i 7 / 12 of Ref. [13]. Interestingly, the observed 
scaling behaviors of density and lengths seem to be in- 
dependent of A as long as A < 1. In next section, we 
introduce our model in detail and present simulation re- 
sults of several values of A. In Sec. Ill, we discuss the 
collective motion of single species domain to understand 
the A-independent kinetics of the reaction. Finally, we 
conclude with summary in Sec. IV. 



II. MODEL AND SIMULATION RESULTS 

On one dimensional lattice of size L with periodic con- 
dition, we consider a bosonic A + B — > reaction with 
mass-dependent fragmentation uniformly driven to the 
right. First, randomly select a site. If the selected site i 
is occupied by mass m, m(= n*) mass hops to the right 
with unit probability. When the opposite species of mass 
rij+i occupies the site i + 1, the reaction rii + i — > rii+i — m 
occurs instantaneously. For m > n.;+i, to — rt; + i mass oc- 
cupies the site The color of mass at the site i + 1 can 
be changed when to > rij+i. If the same species occupies 
the site i + 1, the coagulation of — > Jif+i + to takes 
place. 

We use the following algorithm for determining the 
mass fragmented from a site i. For A > 0, the mass 
fragmented from site i (nf) is not an integer number. In 
that case, integer [nf] mass hops and then unit mass hops 



with probability n\ 



[x] denotes the integer number 



not greater than the real number x. For A < 0, only unit 
mass attempts to hop with the probability n^. With 
random initial distributions of pa(0) = Pb(0) = 1/2, 
we perform Monte Carlo simulations on a chain of size 
L = 1 X 10 7 for A = -1 and L = 3 x 10 6 for other values 
of A up to t = 10 7 . We average densities and various 
lengths over from 50 to 150 independent runs. 

In Fig. [21 we plot the density of A species 
(pA{t)) and its effective exponent defined as — a(t) = 
log[p A {t)/pA(t/b)]/\ogb with b = 5 for A = 
1/2,0,-1/2,-1,-2 respectively. For comparison with 
HC particles, we also plot the results of HC particles. 
With pa{0) = Pb(0), the density Pa($) decays as pa ~ 
i -Q , and decays faster than that of HC particles in all 
stages for A < 1. Especially in early stage, the density 
decay becomes faster as A decreases for A < 0, which 
is reflected by bumps in a(t). After the early stage, 
PA(t) slowly gets into the asymptotic scaling regime of 
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FIG. 2: (Color online) (a) Double logarithmic plot of pA- (b) 
The effective exponent a(t) and the dashed line is a — 1/3. 
The color of each line is same in both panel. 



a = 1/3. From the scaling plot of t a p(t) against t, we es- 
timate a = 0.329(2), 0.331(3), 0.333(2) and 0.331(2) for 
A = 1/2, 0, —1/2 and —1 respectively which agree well 
with previous studies [0, [H, [II [3. For HC particles, 
we estimate a = 0.327(3). Hence our model approaches 
to the scaling regime faster than HC particle system does. 
For A = —2, we estimate a = 0.37(1). It indicates that 
the asymptotic scaling region is not reached yet. The 
physical argument for the slow approach to the asymp- 
totic scaling of the A = — 2 case is discussed in the fol- 
lowing paragraph. 

To understand the early decay, we measure the maxi- 
mum number of particles occupying a single site at time 
t, N max (t). By definition, we have n^t) < N max (t). Fig. 
[3] shows the plot of N max (t) versus t for several A values. 
In early time, like particles pile up on sites due to multiple 
occupation. Hence N max increases to a certain maximum 
value, and after then decreases monotonously. The max- 
imum of N max and the time it takes to reach the value 
tend to diverge as A — * —00. Hence the time it takes to 
reach the asymptotic scaling regime of N max ~ 0(1) is 
also diverge as A — ► —00. For example, the maximum of 



N„ 



for A 



2 is reached at t ~ 3 x 10 . It means that 
■ x ~ 0(1) is very 
Hence the system 



the asymptotic scaling regime of N, 
slowly approached after t > 3 x 10 7 
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FIG. 3: (Color online) Semi-logarithmic plot of N ma x versus 
t. The main plot shows N max of A = —1, —2. In the inset, 
N m ax of A = 1/2, 0, -1/2 are plotted. 



already feels the finite size before t = 10 7 , which results 
in the larger estimate of a. For this reason, we perform 
simulations only for A > — 2. 

Since the drift velocity of a particles on a site with large 
mass is smaller than that on a site with small mass, the 
annihilations of adjacent AB pairs are accelerated. For 
A < 0, it can be easily understood by considering the 
relative drift velocity of adjacent opposite species. When 
the masses of two adjacent sites i and j (j < i) occupied 
by opposite species are rii and rij respectively, the relative 



velocity of a particle at site i to j is Vi 



(1 - (rit/rij) x )/n l . For m > Tij 



is negative so the 



reaction is accelerated. On the other hand, for rii < rij, 
a particle on the site j is hard to react with a particle 



on the site i due to 



> 0. For vn = of the m = 



rij case, diffusive motions lead to the reaction as in HC 
particle system. As a result, the reactions are much more 
enhanced between domain boundaries of t?y < 0, and also 
more accelerated as A — > —00 in early stage. Although 
Nmax decreases as the scaling regime is approached, the 
enhanced reactions by the mass difference should be still 
dominant due to N max > 1. For this reason, the densities 
of A < 1 decay faster than that of HC particles in the 
whole stage. 

Next, we present simulation results of various lengths 
characterizing the spatial organization of masses. The 
average length of a single species domains(^) is defined 
as the distance between the first mass of two adjacent 
opposite species domains. The length Iaa and Iab are 
defined as the inter-mass distance between two adjacent 
masses of the same species and of opposite species re- 
spectively 0]. The lengths algebraically increase in time 
t for p A {0) = p B (0) as 
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FIG. 4: (Color online) The scaling plot of £aa against t. 
The main plot shows the scaling plot of £AAt~ 1 ' ZAA ■ We 
use 1/zaa = 1/3 for A = 3/4, HC particles and 0.332 for 
A = 0, — 1. We add —0.62 to the scaled value of A = —1 and 
—0.02 to that of HC particles for the better presentation. The 
inset shows the double logarithmic plot of Iaa against t. 



FIG. 6: (Color online) The scaling plot of £ against t. The 
main plot shows the scaling plot of lt~ x ' z . We use 1/z — 0.65 
for A = 3/4, -1, HC particles and 0.64 for A = 0. The inset 
shows the effective exponent l/z(t). The solid and dashed 
horizontal lines correspond to 1/z = 2/3 and 7/12 respec- 
tively. 
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FIG. 5: (Color online) The scaling plot of £ab against t. 
The main plot shows the scaling plot of £ab£~ 1 ^ zab ■ We 
use 1/zab = 3/8 for A = 3/4, 0, HC particles and 0.37 for 
A = —1. We add a constant 0.2, 0.4 and 0.5 to the scaled 
value of A = 3/4, 0, and —1 respectively to avoid overlap. 
The inset shows the double logarithmic plot of £ab against t 



In the measurement of lengths, the periodic boundary 
condition is not used for better statistics. Fig. [4] shows 
the scaling plots of l AA t~ x l ZAA for A = 3/4,0,-1. We 
obtain the best scaling plots with 1/zaa = 1/3 for A = 
3/4 and 0.332 for A = 0, —1 respectively. For HC particle 
system, the best scaling plot is obtained with 1/zaa — 
1/3 as expected. For the length Iab (Fig. 0, we obtain 
the best scaling plot of i A B~t~ 1/zAB with 1/z A b = 3/8 
for A = 3/4, and 0.37 for A = — 1 respectively. For 
HC particles, 1/zab = 3/8 gives a nice scaling plot as 
expected. The simulation results of Iaa and £ab agree 



very well with the prediction 1/zaa = 1/3 and 1/zab — 
3/8 of previous studies (l3l. [Til]. 

For the domain length £, 1/z — 2/3 is predicted by 
Ref. [13], while 1/z = 7/12 is predicted and are numeri- 
cally supported by Ref. [HI, [lj] . Our simulation results 
support 1/z = 2/3 rather than 7/12. Fig. [6] shows the 
scaling plots of £t- x l z for A = 3/4,0, -1 and HC parti- 
cles. We obtain the best scaling plots with 1/z = 0.65 for 
A = 3/4, —1, HC particles and 0.64 for A = respectively. 
These values are close to the prediction 1/z = 2/3 of Ref. 
(T3 |. The inset shows the effective exponent 1/ z(t) oi£{t) 
defined as ln(£(bt)/l(t))/\nb with 6 = 5. The effective 
exponent shows that our model reaches asymptotic re- 
gion faster than HC particle system does. On the other 
hand, 1/z = 7/12 indicates the violation of the relation 
p r-j l/\fi in the drift case, while 1/z = 2/3 satisfies the 
relation. We check the relation of p ~ l/Vl by investi- 
gating the scaling plot of p£ x . We obtain the best scaling 
plot with x — 0.51(2) for several A and HC particle sys- 
tem (not shown) , which numerically confirm the relation 
p ~ 1/VI in this uniformly driven case. With a = 1/3, 
we also expect 1/z = 2/3 from the relation p ~ 1/Vi- 
Therefore our numerical results agree well with the pre- 
diction of 1/z = 2/3 [13 • 

Our simulation results indicate that the kinetics of the 
reaction with mass-dependent fragmentation does not 
depend on the exponent A(< 1). In particle version, 
772 = particles hop to the right at the same time in our 
model, while only the boundary particle can hop in HC 
particle system. The fragmented mass m can be regarded 
as TTi particles selected in various ways such as random 
selection or ordered selection from the top of the pile with 
Hi particles. Hence no direct mapping between our model 
and HC system is available. For the ordered BLG model 
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of Ref . [23| , it was shown that the dynamics of the den- 
sity described by Burgers equation and the dispersion of 
the center of mass (CM) scales as i 4 / 3 which is the char- 
acteristic of asymmetric exclusion process (ASEP) [24j . 
For the same kinetics of our model as that of HC system, 
collective motions of single species domain should share 
the feature of ASEP. However, to the best of our knowl- 
edge, there are no such studies on bosonic model with the 
present kind of mass-dependent fragmentation. In next 
section, we discuss the collective diffusion of driven single 
species domain with the mass-dependent fragmentation. 



III. THE COLLECTIVE MOTION OF A SINGLE 
DOMAIN 



s 

u 




For simplicity, we consider the bias to the right only. 
Hence a particle hops to the right with unit probabil- 
ity. With initial density of p, particles are randomly dis- 
tributed on a ring of size L. When randomly selected site 
i is occupied by rii particles, to = nf particles randomly 
selected and moves to the right at the same time. When 
the site i + 1 has rti+i particles, the mass coagulation of 
rit+i — > Uj+x +m takes place instantaneously. For A = 1, 
all rii particles move and only coagulation occurs without 
fragmentation. 

We measure the dispersion of the displacement of the 

center of mass (cr CM ) defined as <j C m = X 2 cm — X CM . 
Xcm is the displacement of the center of mass defined 
as Xcm — X)f X i/N, where N is the total number of 
particles and Xi is the displacement of ith particle from 
its initial position. The dispersion a CM is for driven 
bosonic particles and scales as t 4 ^ 3 for driven HC particles 
[24j . We perform Monte Carlo simulation for various 
values of A from 1/2 to —1 on chains of size L = 10 4 
with periodic boundary condition. We average o CM over 
3 x 10 4 independent runs up to t = 10 4 time steps. In 
simulations, we use the exact p and set p — 0.2 for A = 
1/2, 0.5 for A = and 0.01 for A = — 1 respectively. 

Fig. [7] shows the plot of <r CM t~~ f L against t. <J CM 
crossovers to t 1 scaling regime. We obtain the best scal- 
ing plot with 7 = 1.38(1) for A = 1/2, and 1.37(1) for 
A = — 1 and respectively which is slightly larger than 
7 = 4/3 of HC particles. In the ordered BLG model of 
Ref. [23| , the velocity of CM is given as vcm = 1/(1 + /°) 
for the extreme drift to the one direction. For A = with 
p = 0.5, we estimate vcm — 0-669(1) which also agree 
well with the prediction of Ref. j23|. For A = with 
other densities, we also confirm the prediction. Hence we 
are convinced that the A — case shares the same feature 
of the ordered BLG model in spite of rather large numer- 
ical values of 7. On the other hand, for A ^ 0, there are 
no theoretical predictions for vcm and o C m- However 
the exponent 7 indicates that the collective motions of 
A < 1 close to that of HC particles rather than that of 
bosonic particles. 

We conclude that the collective driven motions of 
masses with mass-dependent fragmentation share the fea- 



FIG. 7: (Color online) The scaling plot of ac M t 1 L against 
t. We use 7 = 1.38 for A = 0.5, and 1.37 for both A = -1 and 
with system size L — 10 4 . We multiply the line of A = — 1 
by 0.04 to draw the line near the others. 



ture of HC particles as long as A < 1. It results in A- 
independent kinetics of A + B — > reaction with the 
mass-dependent fragmentation. 



IV. SUMMARY 

In summary, we investigate the kinetics of uniformly 
driven A + B — > reaction with mass-dependent frag- 
mentation. In this model, the fragmented mass to of a 
site with mass m is given as to = nf, and it is driven 
to the one direction. When opposite species mass occu- 
pies the same site, mass reaction takes place instanta- 
neously. Since the fragmented mass to of A < 1 is less 
than the whole mass rii of a site, the exponent A con- 
trols the strength of the attractive interaction between 
bosonic particles occupying the same site. The A = 
corresponds to the uniformly driven hard-core (HC) par- 
ticles system because only one particle hops with unit 
rate. Compared with A = 0, the attractive interaction of 
A > is weaker than HC interaction, while the interac- 
tion of A < is stronger. Our model generalize the HC 
interaction and enables to investigate how the kinetics of 
the reaction depends on the interaction strength. Also 
we expect that our model gets into the asymptotic scal- 
ing regime faster than HC particle system does, since HC 
constraint is somewhat relaxed in our model by allowing 
any particle on a site to move. 

We perform Monte Carlo simulations in one dimension 
for several values of A, and measure the critical exponents 
of density and various lengths which exhibit power-law 
scaling for pa(0) = Pb(0). The scaling behaviors of den- 
sity and lengths except the domain length £ are the same 
as those of HC particles of Ref. [ill, EE EJ] ■ The scal- 
ing behavior of £ are close to i 2 / 3 of Ref. [l2| rather 
than t 7 / 12 of Ref. [l3j |. Interestingly, the kinetics of the 
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reaction is independent of A as long as A < 1. To under- 
stand A-independent behavior, we investigate the collec- 
tive driven motions of single domains. The displacement 
of the center of mass (a% m) scales as V with 7 = 4/3 
for driven HC particles [24| and <J% M is zero for bosonic 
particles. For various A values, we find that 7 is slightly 
larger than but close to 4/3. The collective motion of a 



single domain shares the feature of HC particles. Fur- 
thermore the motion is independent of A, which results 
in A-independent kinetics of the reaction. However for 
the quantitative understanding of the collective motion, 
another analytical study is desired. 

This work is supported by Grant No. R01-2004-000- 
10148-0 from the Basic Research Program of KOSEF. 
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